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O Abstract 

O ! 

, We study radiative corrections to third generation scalar fermions into gauge bosons 

Z and . We include both SUSY-QCD, QED and full electroweak corrections. It 
is found that the electroweak corrections can be of the same order as the SUSY-QCD 
^ | , corrections and interfere destructively in some region of parameter space. The full 

| one loop correction can reach 10% in some SUGRA scenario, while in general MSSM, 

the one loop correction can reach 20% for large tan (3 and large trilinear soft breaking 
terms A^. 

\D . 

1. Super symmetric theories predict the existence of scalar partners to all known quarks 
and leptons. In Grand unified SUSY models, the third generation of scalar fermions, t, b, r, 
gets a special status; due to the influence of Yukawa- coupling evolution, the light scalar 
fermions of the third generation are expected to be lighter than the scalar fermions of the 
first and second generations. For the same reason, the splitting between the physical masses 
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of the third generation may be large enough to allow the opening of the decay channels like 
: /2 — > fiV and/or / 2 — » where V is a gauge boson and $ is a scalar boson. 
Until now there is no direct evidence for SUSY particles, and under some assumptions on 
^ ■ their decay rates, one can only set lower limits on their masses [1]. It is expected that 
the next generation of e + e~ machines and/or hadron colliders (LHC and Tevatron) could 
establish the first evidence for the existence of SUSY particles if they are not too heavy. 
If SUSY particles would be detected at hadron colliders, their properties can be studied with 
high accuracy at a high-energy linear e + e~ collider [2]. It is thus mandatory to incorporate 
effects beyond leading order into the theoretical predictions, both for production and decay 
rate, in order to match the experimental accuracy. 

In this spirit, radiative corrections to the decays of SUSY particles have been carried out. 
In particular, the QCD corrections to scalar quark decay into quarks plus charginos or 
neutralinos have been studied in [3], while the full one loop analysis has been addressed in 
[4] and found to have important impact on the partial decay widths of scalar fermions. The 
QCD corrections to the decays of heavy scalar quarks into light scalar quarks and Higgs 
bosons are found to be important [5]. 

Obviously, most of the studies concentrated on the production and decay of light states 
t±,bx and fx, while heavier states received less attention [4, 5, 6, 7]. These heavy states can 
be produced both at LHC and/or at the future e + e~ linear colliders. The decay of the heav- 
ier states third generation scalar fermions is more complicated than the light one. One can 
basically have four set of two-body decays: i) Strong decay for stop and sbottom t 2 —> tg, 
i>2 — ¥ bg : if these decay are kinematically open they are the dominant one. ii) decay to 
chargino and neutralino : f 2 — > fXii f2 — > f xt ■ 

If the splitting between light and heavy third generation scalar fermions is large enough 



1 



we may have the following decays: iii) f 2 — > $° = h°,H°,A°, and f 2 — > f' 1 H ± . iv) 

f 2 - f.Z and / 2 -> Z'^. 

In the MSSM, the decay modes / 2 — >• and /j — > f'jW, if open and under some 

assumptions, may be the dominant one. Note also that in several benchmarks scenarios 
for SUSY searches, the bosonic decay of ij and ft; may be the dominant [8]. For exam- 
ple, in SPS5 scenario the dominant bosonic decay have the following branching ratios [8]: 
Br(b! -> W'ti) = 81%, Br(b 2 -> W~t x ) = 64% and Br(t 2 ^ Z%) = 61%. While in SPS1 
scenario, we have: Br(b 2 -> = 34% and £r(t 2 -> ZH X ) = 23%. 

Here, we provide the complete one loop radiative corrections to f 2 — > f\Z° and /, — > f'jW 
including real photon emission [9], and discuss their effects in combination with the SUSY- 
QCD corrections [10]. We show that SUSY-QCD corrections can interfere destructively with 
electroweak one. 

2. The tree-level decay width for qf — > q^V can thus be written as: 



with n(x,y,z) = (x 2 + y 2 + z 2 -2xy-2xz-2yz) 1 / 2 , g z jj 2 = 1 ^{(l(-Q f s 2 v )sme f cos6 f + 
Q/s&r cos 0/ Bin 0/, g z jj 2 = g z jj x and g w jj >2 = ■j^ca&OfWn.Oy, l( is the isospin. 

At one loop level, the decay qf — > q^V receive contributions from vertex diagrams, gauge 
boson V and scalar fermions qf self energies as well as theire mixings (see [9] for details). 
Note that the transitions between gauge bosons and scalar bosons like W ± -H ± , W ± -G ± , 
Z°-A°, Z°-G° are present. Owing to Lorentz invariance, those mixing are proportional to 
Py momentum; then since the vector gauge bosons W and Z are on-shell transverse, those 
transitions vanish. In what follows we will ignore vector-scalar boson mixing. 
We have evaluated the one-loop amplitudes in the 't Hooft-Feynman gauge using FeynArts 
and FormCalc [13, 14]. The one-loop amplitudes are ultraviolet (UV) and infrared (IR) di- 
vergent. The UV singularities are treated by dimensional reduction and are compensated in 
the on-shell renormalization scheme. We have checked explicitly that the results are identi- 
cal in using dimensional reduction and dimensional regularization. The IR singularities are 
regularized with a small fictitious photon mass S and are eliminated by adding to the one 
loop contribution both real-photon and real-gluon emission [9], fi — > fjV'y and fi — > fjVg. 

Recently, there have been several developments in the renormalization of MSSM. Several 
schemes are available [16]. Here, we follow the strategy of [4, 7] by introducing counter- 
terms for the physical parameters, i.e. for masses and mixing angles, and perform field 
renormalization in a way that residues of renormalized propagators can be kept at unity. 
For SM parameters and fields, we will adopt throughout, the on-shell renormalization scheme 
of Refs. [15]. In the on-shell scheme we use the mixing angle sw (resp cw) is defined by 

= 1 — Myy/Mg (resp Cyy = Myy / M|) . Its counter-term is completely fixed by the mass 
counter-terms of W and Z gauge bosons. 

The extra parameters and fields we still have to renormalize in our case are the scalar 
fermion wave functions ft and the mixing angle 0/ which enter in the tree level amplitude 
eq. (1). 
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In the general case, where sfermions mixing is allowed, the wave functions of the two 
sfermions mass eigenstates are not decoupled. Taking into account the mixing, the renor- 
malization of the sfermions wave functions and the mixing angle Of can be performed by 
making the following substitutions in the Lagrangian 

h - Z\{*f x + Zll 2 f 2 , f 2 - Zli 2 f 2 + Zl[ 2 h , f ^9 f + 59 f (2) 

To fix all the renormalization constants, we use the following renormalization conditions: 

i) The on-shell conditions for mw, rnz, m e and the electric charge e are defined as in the 
Standard Model [15]. 

ii) On-shell condition for the scalar fermion /, : we choose to identify the physical scalar 
fermion mass with the corresponding parameter in the renormalized Lagrangian, and require 
the residue of the propagators to have its tree- level value, i.e., 

5Z U = -^{ — (£^( P 2 ))}| p2=m a, , 5Z l3 = J*_ J , 5m\ = ^(S^~(m|)) (3) 

fj fi 

where J^ttIp 2 ), h j — 1>2 is the scalar fermion bare self-energy. 

Jit 3 

iii) For the renormalization condition which defines the mixing angle 9f, we select this condi- 
tion in such a way to kill the transitions fi «-> fj at the one-loop level. The renormalization 
of the scalar fermion mixing angle is then given by [4]: 

60 f = — % ^ — f — (4) 

J 2 m 2 ~ - m 2 ~ 

fj fi 

3. Now we are ready to present our numerical results both for the tree-level and one-loop 
decay widths and branching ratios for fi — > fjZ and fi — > f'jW ± . Let us first fix our inputs 
and SUSY parameters choice. 

As experimental data points, the following input quantities enter: a" 1 = 137.03598, mz = 
91.1875 GeV, mw = 80.45 GeV. For quarks masses, we use effective quark masses that 
reproduce the hadronic vacuum polarization contribution Aa(m|) with a sufficiently high 
accuracy have been chosen [17]. 

For the SUSY parameters, we will use MSSM inputs which look like some of the Snow-mass 
Points and Slopes (SPS) and benchmarks scenarios for SUSY searches [8]. For our study we 
will use SPS1 and SPS5 scenario. As we explained in the introduction, for those 2 scenarios 
the bosonic decays of scalar fermions fi — > fjV, when open, are dominant. 

In SPS1, we have the following spectrum (are listed only the parameters needed): tan (3 = 
10, m A o = 394 GeV, A t = -431.34 GeV, A b = 582.67 GeV, M = 193 GeV, M' = 99 GeV, 
fi = 352 GeV. The mass of the first and second generation scalar fermion is of the order 177 
GeV (average). While the masses of the third generation scalar fermions are : m-jr = 396.43 
GeV, mj 2 = 574.71 GeV , = 491.91 GeV , = 524.59 GeV. The mixing angle are 
given by cos9 t = 0.57 , cos9 b = 0.88. 

In SPS5 (light stop scenario), we have the following spectrum : m^o = 694 GeV, tan j3 = 5 
A t = -785.57 GeV, A b = -139.11 GeV, M = 235 GeV, M' = 121 GeV, fi = 640 GeV The 
mass of the first and second generation scalar fermion is of the order 231 GeV (average). 
The masses of the third generation are m~ tl = 253.66 GeV , m 4 ~ 2 = 644.65 GeV , m bx = 535.86 
GeV , mi 2 = 622.99 GeV. The mixing angle are given by cos 6^ = 0.54 , cos9 b = 0.98. 
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Figure 1: Branching ratios of bosonic decays of b\ (upper plots), b\ (middle plots) and ti 
(lower plots) in SPS1 (left) and SPS5 (right) as function of cos9 t . 



In fact, our strategy is the following : the SPS1 and SPS5 outputs are fixed as above, 
but we will allow a variation of the mixing angles cos 6^, cos 9 b from their SPS values. Ac- 
cording to our parametrization defined in section 2, we choose as independent parameters 
mj 2 ,m^m^, 9 t , 9 b together with /i and tan (3. and Af are fixed by 

m~ = 7T7r( m w cos 23 — m~ sin 2 9 t + mf- sin 2 9 b + m~ cos 2 9 b + mf — m?) (5) 

fl cos 2 ^ * 2 62 bl 

2 2 
m~ — m~ 

A f = [i(tan(J)- 2I f + — ^sin^cos^ (6) 

The variation of cos#( and cos#& imply the variation of mj 2 as well as A t and A b . Since we 
allow variation of the cos^t and mj 2 mass, our inputs can be viewed as a general MSSM 
inputs and not as SPS one. 

As outlined in section 2, A t ^ are fixed by tree level relation eq. (6). Of course, A t> b receive 
radiative corrections at high order. However, A t and A b enter game only at one-loop level 
in our processes, radiative corrections to A t and A b is considered as two-loop effects. 
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Figure 2: Relative correction (electroweak EW, SUSY-QCD and full) to b 2 — > t{W as func- 
tion of co$0 t in SPS1 (left) and SPS5 (right) 
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Figure 3: Relative correction to b\ — > t±W and t 2 — > t\Z as function of cos# t in SPS1 (left) 
and SPS5 (right) 



In Fig. (1) we show branching ratios of b±, b 2 and t 2 . We evaluate the bosonic decays : 
bi — > W~ti, b 2 — > W~t\ and t 2 — > Z°ti as well as the fermionic decays fa — > Xifixtf) 
as function of cos9 t for SPS1 (left) and SPS5 (right) scenario. From those plots, it is clear 
that the bosonic decay, once open, are the dominant one for |cos#t| ~ 0.4 — > 0.45. For 
| cos#j| ~ 0.4 the light stop vnq is about 100 GeV, when | cos^| increases, the vnq increases 
also and for large | cos 0* | the bosonic decays are already close and the branching ratio 
vanishes. 

We note that in the case of SPS1 the bosonic decays are open only for 0.4 < | cos 6^1 < 0.6 
Fig .(1) (left). In the region | cos 6^1 < 0.4, the light stop is below the experimental upper 
limit ~ 90 GeV, and no data are shown. While in the case of SPS5 Fig .(1) (right), for 
| cos^| < 0.5, we find that mj is below the experimental upper limit and also 6p > 0.001 
due to large splitting between stops and sbottoms. 

The magnitude of SUSY radiative corrections can be described by the relative correction 
which we define as: 

A = r^(f t -+ JjV) - t*™{1 -> JjV) 
r tree (/, ^ fjV) 

In Fig. (2) we illustrate the relative correction A as function of cos 6^ for the decay 
b 2 — > Wt\ in SPS1 (left) and SPS5 (right). As it can be seen from the left plot, the SUSY- 
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QCD corrections [18] lies in the range — 1% — * —6% while the EW corrections lie in the range 
4% -> 10% for cos0 t « 0.4 -> 0.65. The SUSY-QCD and EW corrections are of opposite 
sign, there is a destructive interference and so the full one-loop corrections lie between them. 
For cos9 t — > 0.65, the total correction increases to about 10%. This is due to the fact that 
for cosOt — > 0.65 the mass of light stop is mj ~ 444 GeV, the decay b 2 - ► Wt\ is closed 
and so the tree level width decreases to zero. The observed peaks around cos9 t ~ 0.46 
(resp cos^t ~ 0.53) correspond to the opening of the transition ti — > Xit (resp t\ — > xfy)- 
The right plot of Fig. (2) in SPS5 scenario, exhibits almost the same behavior as the left 
plot. The electroweak corrections interfere destructively with the SUSY-QCD ones, the full 
corrections are between —4% — > 10% for cos #4 G [0.5,0.9]. 

In Fig. (3) we show the relative correction A as function of cos# t for the decay b\ —>■ Wt\ 
and t 2 Zti in SPS1 (left) and SPS5 (right) scenario. 

In the case of t 2 — ► Zt\, the total correction lies in — 1 — > 7% (resp —5 — > 6%) in SPS1 
(resp SPS5) scenario. From Fig. (3), one can see that the relative corrections for b± — > Wt\ 
are enhanced for cos #4 ~ 0.6 (resp cos#t ~ 0.75) in SPS1 (resp SPS5). This behavior has 
the same explanation as for b 2 — * Wt\ in figure. (2). At cos# t ~ 0.6 (resp cos# t ~ 0.75) in 
SPS1 (resp SPS5), the decay channel b\ — > Wt\ (resp t 2 — > Zt\) is closed and so the tree 
level width decreases to zero. The observed peaks around cos# 4 ~ 0.46 (resp cos 6^ ~ 0.53) 
correspond to the opening of the transition t\ — > xl^ ( r esp t\ — > 

In all cases, we have isolated the QED corrections (virtual photons and real photons), we 
have checked that this contribution is very small, less than about 1%. 

Fig. (4) illustrates the relative corrections to t 2 — > b\W, t 2 — > t\Z (left) and b 2 — > b±Z, 
b 2 — > t\W (right) as function of Ab = A t in general MSSM for large tan (3 = 60, \i = 500 GeV, 
M 2 = 130 GeV and Ma = 200 GeV. It is clear from this plot that the relative corrections 
are bigger than in the cases of SPS scenarios. This enhancement shows up for large \Ab\ 
and also near threshold regions. In this scenario, the SUSY-QCD corrections are about 
2%, the electroweak corrections are about 5% while the QED corrections are very small. 
The dominant contribution comes from the Yukawa corrections and is enhanced by large 
tan/3 = 60 and large \Ab\. 

In the left plot of Fig. (4), the region \Ab\ = \A t \ < 300 GeV has no data. This is due to 
the fact that splitting between t 2 and t\ (t 2 and bi) is not large enough to allow the decays 
t 2 — > t\Z and t 2 —> b\W . 

In the right plot of Fig. (4), when \A b \ = \A t \ w GeV, the splitting between b 2 and t\ is 
close to mw mass and so the tree level width for b 2 — > t{W + almost vanish, consequently 
the correction is getting bigger. This behavior has been also observed in previous plots for 
b 2 -> W+. 

Finally, in Fig. (5) we illustrate the decay width of t 2 — > t\Z as function of mq in SPS1 
(left) and SPS5 (right). In SPS1 (resp SPS5) the decay width of t 2 -> t x Z is about 8 GeV 
(resp 13 GeV) for light stop mass of the order 100 GeV. Obviously, these decays width 
decrease as the light stop mass increase. 

It is clear that the SUSY-QCD corrections reduces the width while the electroweak cor- 
rections cancel part of those QCD corrections. Both in SPS1 and SPS5, the full one loop 
width of t 2 —>t\Z is in some case slightly bigger than the tree level width. 

4. To conclude, a full one-loop calculations of third-generation scalar-fermion decays into 
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Figure 4: Relative correction to t 2 — » b\W,t 2 — > t\Z (left) and b 2 — > b\Z ,b 2 —> t\W (right) 
as function of A t = A b in MSSM for /i, M 2 , M A = 500, 130, 200 GeV and tan (3 = 60 




Figure 5: Tree and one loop decay width of t 2 — > t x Z as function of 

gauge bosons W and Z are presented in the on-shell scheme. We include both electroweak, 
QED and SUSY-QCD contributions to the decay width. It is found that the electroweak 
and SUSY-QCD corrections interfere destructively. 

The size of the one-loop effects are typically of the order —5% — > 10 % in SPS scenarios which 
are based on SUGRA assumptions. While in the general MSSM, the size of the corrections 
are bigger and can reach about 20% for large tan/3 and large soft SUSY breaking At,. Their 
inclusion in phenomenological studies and analyses are then well motivated. 
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